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RADIATIVE  PREHEAT  IN  STRONGLY  COUPLED, 
LASER  ACCELERATED  PLASMAS 


Introduction 

Control  of  the  Rayleigh-Taylor  instability  in  ablatively  driven  laser  targets  is  one  of  the 
primary  challenges  in  achieving  high  gain,  direct  drive  laser  fusion.  One  potential  mechanism 
for  inhibiting  the  instability  is  the  deposition  of  high  energy  plasma  radiation  at  die  ablation 
front  in  order  to  heat  and  thereby  alter  the  density  profile.  A  moderate  atomic  number  material 
doped  into  the  plastic  target  will,  when  heated  in  the  corona,  produce  the  penetrating  x-rays 
absorbed  by  the  dense  layer.  Gardner,  et  al.[l]  suggest  that  the  ablative  Rayleigh-Taylor  formula 
for  the  growth  rate  7  could  be  written  as 


(1) 


where  k  is  die  unstable  wavenumber,  g  is  the  effective  acceleration,  pabi  is  Ac  density  at  the 
ablation  surface,  A  =|  Vp/p  )««,  and  m  is  the  ablative  mass  flux.  Stabilization  might  be 
achieved  if  the  density  as  well  as  its  gradient  can  be  reduced  at  the  ablation  front  by  radiative 
preheat. 

The  choice  of  which  dopant  and  how  much  to  add  to  the  target  is  central  to  this  proposed 
solution  for  damping  hydrodynamic  instabilities.  Due  to  the  large  range  of  possibilities, 
numerical  simulations  offer  a  viable  means  of  discrimination.  However,  Acre  are  important 
physical  issues  which  need  to  be  considered  in  order  to  make  a  proper  assessment,  (a)  Due  to 
Ae  high  sound  speed  and  electron  degeneracy  pressure  of  solids  a  realistic  equation-of-state 
(EOS)  is  required,  (b)  The  large  range  in  density  and  temperature  from  Ae  emitting  corona 
to  Ae  absorbing  target  suggest  Ae  use  of  non-LTE  ionization  dynamics  and  radiative  transfer 
rather  than  diffusion,  (c)  In  high  density  regions  such  as  Ae  compressed,  accelerating  target,  Ae 
plasma  can  become  strongly  coupled.  The  importance  of  this  last  effect  is  typically  estimated 
through  Ae  plasma  coupling  parameter: 

r- JlV-  (  n'  V/8  roi 

r  kBTRJS  ~5revU022cm~3i  ’ 

where  e  is  Ae  electron  charge,  kB  is  Boltzmann’s  constant,  T  is  Ae  temperature,  Ris  is  Ae 
ion  sphere  radius  for  Ae  ion  density  n<,  and  Z  is  Ae  mean  charge  state.  For  T  >  1  Ae  plasma 
ionization,  radiation,  and  opacity  can  be  altered,  which  in  turn  will  alter  Ae  plasma  dynamics. 
Furthermore,  Ae  spectrum  of  Ae  radiation  penetrating  to  Ae  fuel  may  be  affected  and  this  is 
relevant  to  unwanted  fuel  preheat. 

As  part  of  a  program  to  assess  Ae  radiative  preheat  effects  NRL  has  been  developing  a 
1-D  Lagrangian  radiation  hydrodynamics  code  which  includes  the  above  physical  processes  for 
Manuscript  approved  March  2,  1993. 


a  range  of  materials.  In  the  present  report  we  will  compare  the  density  evolution  of  the  target  as 
the  complexity  of  the  model  employed  is  increased  in  three  stages:  (a)  a  high  density  EOS,  (b) 
plus  a  collisional-radiative  equilibrium  (CRE)  ionization  dynamics  and  probabilistic  radiative 
transfer,  (c)  plus  strongly  coupled  plasma  effects.  In  each  case  the  planar  target  is  60  pm  thick 
and  composed  of  50%  carbon  and  50%  hydrogen.  The  incident  laser  corresponds  to  the  KrF 
NIKE  laser  system  with  0.264  p m  light  and  an  early  version  of  the  incident  pulse  as  displayed 
in  Figure  1.  The  particular  choice  of  the  dopant  material  and  model  improvements  will  be 
addressed  in  future  work. 


Case  (a):  High  Density  EOS 

Let  us  first  consider  the  target  evolution  when  a  high  density  EOS  is  coupled  with  the 
hydrodynamics,  but  radiative  losses  are  neglected.  The  EOS  consists  of  a  Cowan  ion  model 
using  the  Debye  and  melt  temperatures,  the  Thomas-Fermi  electron  cloud  with  Z -scaling  for 
arbitrary  compositions,  and  a  Barnes  correction  to  obtain  the  proper  pressure  during  initial 
compression  from  the  cold  solid.[2]  The  thermal  conductivity  for  the  dense  plasma  includes 
electron  degeneracy  and  corrections  to  the  Coulomb  logarithm  according  to  Lee  and  More. [3] 
Figure  2  shows  the  evolution  of  the  density  profile  with  the  characteristic  highly 
compressed,  rearward  accelerating  region  behind  the  ablation  front.  The  steep  density 
profile  leading  to  the  ablation  front  at  peak  density  is  highly  Rayleigh-Taylor  unstable. 


Case  (b):  Ionization-Dynamics  and  Radiation 

The  CRE  model  for  the  atomic  populations  uses  rates  for  collisional  and  photo-ionization, 
3-body,  dielectronic,  and  radiative  recombination,  collisional  excitation  and  photo-pumping, 
collisional  de-excitation  and  radiative  decay,  and  inner  shell  photo-absorption.  [4]  The  detailed 
configuration  accounting  includes  the  stripped  state  for  C  and  H,  the  ten  principal  quantum 
levels  for  H,  six  ground  states  and  48  excited  levels  for  C.  A  probability-of-escape  formalism  is 
used  for  the  bound-bound  and  bound-free  radiation  transport, [5]  and  a  multi-frequency  approach 
for  the  free-free  radiation.  Detailed  comparisons  of  the  EOS  and  CRE  calculations  established 
a  demarcation  line  in  the  p-T  plane  across  which  a  smooth  transition  form  the  EOS  charge  state 
and  pressure  to  the  CRE  values  can  be  maintained. 

Hie  evolution  of  the  density  profile  in  Figure  3  shows  in  comparison  with  case  (a)  a 
significant  reduction  in  p  and  |  Vp  |  at  the  ablation  front,  and  an  increase  in  the  separation 
between  the  critical  surface  and  the  ablation  front.  The  critical  surface  occurs  roughly  at 
the  density  0.3  gm/cms.  In  the  present  case,  |  Vp/p  j abt  is  reduced  between  5  and  10,  and 
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by  eqn.  (1),  the  critical  wavelength  for  instability  would  move  to  larger  values.  However, 
calculating  T  from  eqn.  (2)  at  the  ablation  surface  shows  that  it  is  over  10  at  the  peak  of  the  laser 
pulse,  and  remains  above  2  even  as  late  as  8  nsec.  This  suggests  that  strong  coupling  effects 
should  be  taken  into  account  and  we  turn  to  this  model  next. 
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Case  (c):  Strong  Coupling  in  Laser  Plasmas 

High  density  effects  on  plasma  ionization  dynamics  has  been  an  active  area  of  research 
for  several  decades,  yet  a  practical  theory  applicable  in  hydrodynamic  simulations  is  still  under 
development.  The  model  we  employ  is  based  upon  the  work  of  Hummer  and  Mihalas.[6] 
The  decrease  in  the  population  of  a  particular  atomic  level  is  accounted  for  smoothly  and  in  a 
thermodynamically  consistent  manner  through  a  reduction  in  the  effective  statistical  weight  gni 
'  of  that  level.  Pressure  ionization  is  treated  as  an  excluded  volume  correction  through 

X 

9ni  =  S'niCxp{-Jr[rM.i(n)/JR/5]8}, 


where  rorb( n)  is  the  bound  electron  orbital  radius  for  the  level  with  principal  quantum  number  n 
and  f  is  the  fraction  of  neutral  atoms.  The  microfield  perturbations  to  the  nl  level  of  an  ion  of 
charge  state  Z»  is  assumed  to  follow  a  Holtsmark  distribution  resulting  in  the  further  reduction 


Cl  -  <7n/Gxp{-16( 


y/(Zj  +  1  y 

Xnl 


where  Xnl  is  the  ionization  potential  of  the  level.  To  account  for  the  electron  screening  in 
strongly  coupled  plasmas  we  include  a  continuum  energy  lowering  term  which  transitions  from 
the  Debye-Huckel  to  die  ion-sphere  forms: 


A  E,ov  - 


_ Zf _ 

max{2Rrs/3>RDHZ/(Zi  + 1)* 


where  Rdh  is  the  Debye-Huckel  radius. 

Besides  changing  the  effective  statistical  weight  of  a  level  one  must  also  reduce  the  atomic 
physics  rates  to  ensure  that  the  proper  equilibrium  is  attained  at  high  densities.  We  follow  the 
procedure  of  Busquet[7J. 

The  density  evolution  for  the  full  model  is  shown  in  Fig.4.  The  most  surprising  feature  is 
the  deterioration  of  the  dense  slab  near  the  rear  of  the  accelerating  target  and  the  consequent 
broad  density  profile.  The  strong  coupling  has  reduced  the  allowed  population  in  so  many  levels 
that  the  internal  energy  can  no  longer  act  as  a  heat  sink.  As  a  result,  the  rearside  temperature 
and  charge  state  are  highest  in  case  (c). 
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Another  interesting  feature  is  the  radiation  spectrum  passing  through  the  rear  side  of  the 
target  near  peak  laser  intensity  for  cases  (b)  and  (c)  shown  in  Fig .5.  Note  the  stairstep  feature  of 
the  emitted  bound-free  continuum  radiation  in  case  (b)  reduces  to  fewer  stairs  when  the  strong 
coupling  effects  are  included.  Furthermore,  the  mean  photon  energy  is  smaller  for  case  (b). 
This  difference  in  the  spectrum  is  significant  for  heating  of  the  fuel.  The  fraction  of  the  DT 
fuel  which  is  neutral  will  be  most  opaque  to  ~  IS  eV  photons  and  the  absorbtion  cross-section 
decreases  as  the  energy  cubed  above  the  threshold.  In  the  present  model,  the  inclusion  of  strong 
coupling  effects  drastically  reduces  the  rearside  spectrum  in  the  10  to  100  eV  range  and  implies 
a  reduction  in  the  photo-absorption  heating  of  the  fuel. 

Discussion 

Comparison  of  Figs. 2,  3,  and  4  demonstrate  substantial  changes  in  the  target  density  as 
the  complexity  of  the  simulation  is  progressively  increased.  Case  (b)  indicates  the  importance 
of  radiation  even  for  targets  composed  of  low  atomic  number  material  such  as  carbon.  This 
conclusion  agrees  with  Duston,  et  al.[ 8]  who  studied  the  problem  with  lower  laser  intensities 
and  longer  wavelengths.  The  role  intended  for  doping  by  moderate  Z -material  appears  already 
to  be  fulfilled  by  carbon  in  typical  plastics.  In  case  (c)  the  high  density  effects  of  pressure 
ionization,  microfield  perturbations,  and  continuum  lowering,  on  atomic  processes  further 
accentuates  the  decrease  in  pabi  and  the  broadening  of  the  density  profile.  The  present  results 
of  including  detailed  non-LTE  radiation  and  strong  coupling  effects  suggest  an  inhibition  of  the 
ablative  Rayleigh-Taylor  instability  in  the  accelerating  target  according  to  eqn.  (1). 

Furthermore,  reduction  of  ions  in  excited  states  due  to  high  densities  shifts  the  backside 
spectrum  toward  higher  photon  energies  and  away  from  the  strong  absorption  peak  of  neutral 
DT  fuel.  However,  these  potentially  positive  conclusions  are  offset  by  the  extreme  deterioration 
from  a  dense,  cold,  accelerated  target  [case  (a)]  to  a  more  diffuse  and  hot  one  [case  (c)]. 

Future  work  will  need  to  concentrate  on  improving  the  microfield  distribution  to  account 
for  large  F’s  and  allow  for  frequency  shifts  due  to  continuum  lowering  toward  the  red  end 
in  the  bound-free  emission.  This  last  aspect,  only  feasible  within  a  complete  multi-frequency 
radiation  transport,  may  alter  some  of  the  present  conclusions  on  the  density  profile  and  fuel 
preheat. 
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Fig.l  Profile  of  incident  laser  pulse  used  in  the  modeling. 
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Fig.3  Evolution  of  the  mass  density  at  the  listed  times  for  case  (b)  with  an  EOS,  CRE 
ionization  dynamics,  and  radiation  transfer. 


8 


